Comment on "Optical Response of Strongly Coupled Nanopraticles in Dimer Arrays 

(Phys. Rev. B 71(4), 045404, 2005). 
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I have re-calculated the extinction spectra of aggregates of two silver nanospheres shown in Figs. 2 
and 3 of Ref. 8. I have used the approximate method of images according to Ref. 8 and an exact 
numerical technique. I have found that the three sets of data (those I have obtained by the method of 
images, the numerical results, and the results published in Ref. 8) do not coincide. In this Comment, 
I discuss the reasons for these discrepancies and the general applicability of the method of images 
to the quasi-static electromagnetic problem of two interacting nanospheres. 
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The method of images (MOD is a powerful tool for 



solving electrostatic problems yj. In particular, it has 
been used to calculate the electrostatic force between 
two charged spheres @- It was shown that the force 
can be well approximated by the Coulomb formula when 
the spheres are far apart. However, as the spheres ap- 
proach each other, they can not be effectively replaced 
by point charges and the Coulomb formula overestimates 
the actual force. The MOI was recently extended beyond 
the electrostatics 0, 0, IE IE 

However, this gen- 
eralization was subject of some controversy [9l Il0|. In 
particular, I have argued that the MOI, as used in the 
above references, is not a physically justifiable approx- 
imation and, therefore, can not be used for calculating 
electromagnetic responses of interacting spheres at finite 
frequencies I3. I have further argued that the formulas 
derived in 0, 0, IE IS do not provide accurate results 
even within electrostatics, when the MOI is, in princi- 
ple, applicable. The authors of Ref. argued that the 
MOI is an accurate approximation at sufficiently low fre- 
quencies, e.g., for / < 1GHz (this claim will be discussed 
below). However, in Ref. [8|, which is the subject of this 
Comment, the MOI is used for the Drude dielectric func- 
tion in the spectral range from 250nm to 1500nm, i.e., 
at much larger electromagnetic frequencies. I have recal- 
culated the data shown in Figs. 2,3 of Ref. using the 
MOI as it is defined in H 0, 0, Q, . I have also cal- 
culated the relevant quantities using the exact method 
(e.g., see 0). I have found that the three sets of data 
(i.e., the data shown in Figs. 2,3 of Ref. 0, the data I've 
obtained according the MOI using the same formulas as 
in Ref. 0|, and the exact results) do not coincide. The 
deviations are significant. This and other points rele- 
vant to the applicability of the MOI are discussed in this 
Comment. 

First, we specify the dielectric function used in calcu- 
lations. The expression given in Ref. |8j is 



eh 



(1) 



where 



1.72 x 10 lb rad/sec, e h = 5.45, e s = 6.18, and 



the relaxation constant 7 is size-dependent: 



7 = vp 1 1 + 2vf/cI , 



where vf = 1-38 x 10 6 m/sec is the Fermi velocity, 
i = 52nm is the electron free path and d is the sphere 
diameter. All quantities give good approximation for sil- 
ver. Since the electromagnetic frequency is expressed in 
the units of energy in Figs. 2,3 of Ref. 0, we re- write the 
expression as 



E{E + iT) 



where E = Two, E„ = 9.68eV and 



T = 0.00181 x E P {1 + 2£/d) 



(3) 



(4) 
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Fig. 1. Size-dependent dielectric function for different 
values of the sphere radius d in the spectral range of interest. 

The dielectric function defined by these equations is 
shown in Fig. 1 for different values of the sphere diame- 
ter d. It should be noted that an analogous graph is also 
shown in Fig. 1 of Ref. 8] for d — lOnm. While the real 
part of the dielectric function is qualitatively similar in 
both figures, the imaginary parts (for d = lOnm) are very 
different. By comparison with other curves in Fig. 1,1 
infer that the dashed curve shown in Fig. 1 of Ref. Q 
corresponds to the case d — * 00 rather than to d = lOnm, 
as claimed. 

The mathematical formalism of MOI 0, S S H, S 
can be summarized as follows. The polarizability of each 
sphere in a two-sphere aggregate (radius of each sphere 



2 



is d and the surface-to-surface separation is a) is given 
by 



(d/2f 



71=1 S 



F, 



S L > T ) 



(5) 



where s — £ m /(e - e m ) is the spectral parameter of the 
theory, with e m being the dielectric function of the trans- 
parent matrix in which the spheres are embedded, the 
upper index L corresponds to polarization of the exter- 
nal field parallel to the axis of symmetry of two spheres 
and the index T corresponds to orthogonal polarization. 
The above formula is quite general and is not a subject 
of controversy. The values of F n and Sn can be, in 
principle, found numerically. The potential advantage of 
the theoretical developments of Refs. 0, 0, IS IE is 
that it provides approximate analytical expressions for 
these quantities: 



F n = An(n + 1) sinh 3 aexp[-(2n + I) a] , 
4 L ) = i{l-2exp[-(2n+l)a]} , 

8<p = i{l + exp[-(2n + l)a]} , 



a = In 



a lot a\ 

i+ W^( 2+ rf) 



(6) 
(7) 

(8) 
(9) 



However, I have previously argued that (i) the values of 
these coefficients can not be, in principle, found from the 
MOI, even approximately, and (ii) the above expressions 
are inconsistent with the electrostatic limit There- 
fore, the formulas 10-1® do not provide a physically 
meaningful approximation. This is illustrated in Figs. 2 
and 3 below. 

We note that a slight change of notations has been 
adopted. Thus, the factors F n defined in © differ from 
those of Refs. 0, S H @, H H by the overall factor of -3 
and the spectral parameter s = e m /(e — e rn ) by the factor 
of —1. This corresponds to the more conventional nota- 
tions 

El USUI- In particular, oscillator strengths © 
satisfy the sum rule ^2 n>0 F n = 1. Also, the expression 
for the polarizability given in Ref. || contains an extra 
factor of e m compared to Eq. JSJ. Although the inclu- 
sion of this factor in the expression for the polarizability 
is incorrect, as can be easily seen in the limit a/d — > oo, 
it does not change any of the spectral lineshapes. 

In what follows, we consider only the results for the 
polarization of the external field being parallel to the 
axis of symmetry of two spheres, since the multipolc 
interaction is strongest in this case. I have used the 
dielectric function defined above to calculate the extinc- 
tion cross-section of the bisphere aggregate for the same 
sets of parameters as in Figs. 2,3 of Ref. @. Namely, 
the dielectric constant of the matrix was e m = (1.61) 2 , 
the sphere diameters were chosen to be d — 5nm (Fig. 2) 
and d = lOnm (Fig. 3), and the ratio a/d was 0.1 and 



0.3 (Fig. 2) and 0.05, 0.15, 0.25 and 0.35 (Fig. 3). The 
extinction in Ref. Q was plotted in arbitrary units and 
not defined precisely. Therefore, I plot the quantity 
(E/3)lm^2 n>0 f n /(s + s n ), where E is expressed in 
electron-volts. This quantity differs from the actual 
extinction cross section only by a constant factor, 
and I have found that it has approximately the same 
numerical values as the data shown in Figs. 2,3 of Ref. . 
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Fig. 2. Extinction spectra for a two-sphere aggregate 
in a transparent matrix with e m = 1.61 2 and d — 5nm. 
Polarization of the external field is parallel to the axis of 
symmetry. The exact result are compared to MOI and to the 
spectra for isolated (non-interacting) spheres. 



We first discuss Fig. 2 which is analogous to Fig. 2 
of Ref. Q , except that data for orthogonal polarization 
are not shown. The exact spectra were calculated by 
the method of Ref. [Tlj with the maximum multipole 
order L = 400 and convergence was verified by doubling 
this number. The MOI calculations were carried out by 
truncating summation in © at n — 400, which is much 
more than is needed for convergence. First, it should be 
noted that the shape of the curves that I have obtained 
by the MOI are markedly different from those shown in 
Fig. 2 of Ref. @- In particular, the second, low-energy 
peak in 
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Fig. 3. Same as in Fig. 2, but for d — lOnm and a different 
selections of the ratio a/d. 



the extinction spectrum for a/d = 0.1 is significantly less 



pronounced in Ref. |8| than in my data. A noticeable 
deviation is also visible in the case a/d = 0.3. A pos- 
sible cause of this discrepancy is discussed below. More 
importantly, the MOI curves in both cases differ from 
the exact result. The difference is very apparent at the 
smaller separation {a/d — 0.1), and still visible at the 
relatively large separation a/d = 0.3. 

The spectra in Fig. 2 are characterized by very strong 
relaxation, because the ratio 2£/d is in this case of the 
order of 20. Thus, the finite size contribution to the 
relaxation constant is approximately 20 times larger 
than the respective constant in bulk. We then consider 
Fig. 3 (d = lOnm), where the relaxation is not as strong. 
The results are shown in Fig. 3 which corresponds to 
Fig. 3 of Ref. || with the exception that the results for 
orthogonal polarization are not shown. Again, there 
is a clearly visible difference between the MOI results 
obtained here and in Ref. :S]. In all cases, the MOI 
spectra are very different from the exact spectra. This 
is especially apparent at the relatively small separation 
a/d = 0.05 when the MOI predicts a spectral peak at 
E w I.5Ev which is not present in the exact data. Even 
for the relatively large separation a/d = 0.35, the MOI 
produces a two-peak structure, while the exact spectrum 
has only one peak. (We note here that in Fig. 3d of 
Ref. Q , the respective curve has also only one peak, but 
its maximum is about 10% smaller than the maximum 
of the spectrum in the noninteracting case. In the exact 
result, the maximum is approximately equal to that for 
the non-interacting case.) 
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Fig. 4. The MOI result for the extinction spectrum of two 
spheres of different diameters d and the ratio a/d = 0.35. 
Polarization of the external field is parallel to the axis of 
symmetry. 

The conclusion that can be made so far is that the 
MOI is inadequate for the spectral range and set of pa- 
rameters used in Figs. 2,3 of Ref. The inaccuracy of 
the MOI is especially evident at smaller inter-sphere sep- 
arations and for larger sphere diameters. We now discuss 
the possible cause of the discrepancy of the MOI calcula- 
tions presented here and in Ref. In Fig. 4, we plot the 
MOI curve for a/d = 0.35 and different values of d. The 
two-peak spectrum obtained at d — lOnm is the same 
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as the one shown in Fig. 3d, while the single peak spec- 
trum obtained at d = 5nm closely resembles the curve 
shown in Fig. 3d of Ref. 0. Thus, the possible cause of 
the discrepancy is that in Ref. @ the actual value of the 
sphere diameter used in calculations was twice smaller 
than what is shown in figure captions. That is, calcu- 
lations in Fig. 2 of Ref. Q were actually performed for 
d = 2.5nm and in Fig. 3 of Ref. || for d = 5nm. Under 
these circumstances, the relaxation due to the finite size 
effects is extremely strong and the spectral parameter s 
has a large imaginary part which effectively weakens the 
multipole interaction of the spheres. 

We now discuss the argument given in Ref. |10j that the 
MOI is accurate as long as the denominators in Eq.[5]do 
not vanish and the oscillator strengths F n 's obey the ap- 
propriate sum rule. In general, this statement is correct 
when |s + s„| 1 for all n. But the result obtained from 
the MOI in this limit is, essentally, the non-interacting 
result. Therefore, application of the MOI under these 
conditions is simply not necessary. This is discussed in 
more detail below. 

First, we note that s is a complex variable while the 
depolarization factors s n are all real. The exact val- 
ues of s n satisfy the inequality < s n < 1, while the 

formulas 0,© can result in negative values of if 
a/d < (2 1 / 3 - l) 2 /2 4 / 3 w 0.0268. It can be seen that the 

smallest possible value of that can be obtained from 
formula Q is —1/3. If the complex value of s is suffi- 
ciently separated from all s n 's in the complex plane, one 
can replace the denominator s — s n by s. (In the mean- 
field approximation, the denominators are replaced by 
s — Q, where Q is the appropriate average of the interac- 
tion operator |15|.) The result is (taking into account the 
sum rule for F n 's) the polarizability of an isolated (non- 
interacting) sphere. We can further expand the result in 
powers of the small parameter s n /s and thus obtain cor- 
rections to the non- interacting result. Unlike the former, 
these corrections depend on the particular choice of F n 
and s n . An important point is that even if the corrections 
are small, they are not necessarily physically meaningful. 
I have previously demonstrated Q that the corrections 



to the non-interacting polarizability obtained with the 
particular choice of F n and s n @-|JSJ| are inaccurate for 
the ratio a/d smaller than m 0.3. Thus, any coincidences 
between the MOI and the exact spectra, such as the ones 
shown in Ref. ^(j; are due to the fact that the multi- 
pole interaction is very weak for the particular choice of 
parameters, and the obtained MOI spectra are, essen- 
tially, the spectra of non- interacting spheres. But in the 
spectral regions where interaction is essential, there is 
an obvious discrepancy between the MOI and the exact 
results [Io|. 

We note that the complex spectral parameter s can 
be removed from the section of the real axis occupied by 
the factors s„, in particular, due to strong absorption. 
In that case, s acquires a large imaginary part. This 
was the case for the simulations shown in Figs. 2,3 of 
Ref. || . Here the resonant interaction of the spheres was 
suppressed by finite size effects ©• In addition, the phys- 
ical size of the spheres used in calculations appears to be 
smaller by the factor of 2 than what is claimed in figure 
captions. This resulted in relatively modest changes in 
the spectra (compared to the spectra of non-interacting 
spheres) which appear to be realistic. However, com- 
parison with numerical results clearly demonstrates that 
even if the very strong relaxation is taken into account, 
the MOI results are qualitatively inaccurate. 

Finally, it was suggested in Ref. 10] that at small elec- 
tromagnetic frequencies (e.g., / < 1GHz), the real part 
of the spectral parameter s is always positive, so that 
the denominators s + s n can never vanish. This is clearly 
incorrect for both conductors and dielectrics. In the case 
of conductors, the real part of the dielectric function 
is negative and of large magnitude. Then the real part of 
the spectral parameter e m /(e — e m ) is also negative and 
approaches zero from the left. It should be emphasized 
that in all realistic cases eh is of the order of unity and 
can not compensate for the term —uj^/u} 2 in the spec- 
tral range / < 1GHz. In the case of dielectrics, the 
low-frequency limit of the real part of s can be either 
positive or negative, depending on the sign of e — e m . 
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